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We present and discuss the analytical r-mode solution to the linearized hydrodynamic equations 
of a slowly rotating, Newtonian, barotropic, non-magnetized, perfect-fluid star in which the gravi¬ 
tational radiation reaction force is present. 


I. INTRODUCTION 

The r-modes, which are pulsation modes of rotating 
stars that have the Coriolis force as their restoring force, 
are driven unstable by gravitational radiation (GR) Q. 
In the frame co-rotating with the star, the energy of the 
unstable r-mode grows exponentially, E = 
with the gravitational timescale tgr given by |2, Q 


II. HYDRODYNAMIC EQUATIONS WITH GR 
REACTION FORCE 

The Newtonian hydrodynamic equations for a uni¬ 
formly rotating, barotropic, non-magnetized, perfect- 
fluid star in the presence of the gravitational radiation 
(GR) reaction force are the Euler, continuity, and Pois¬ 
son equations given, respectively, by 


''"GR = - 


52 X 38 
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where G is the Newton’s constant, c is the velocity of 
light, Q is the angular velocity of the star and J = 
drpr^ with p and R being, respectively, the density 
and the surface’s radius of the unperturbed star. For a 
wide range of relevant temperatures and angular veloci¬ 
ties of newly born, hot, rapidly-rotating stars, bulk and 
shear viscosity do not suppress the exponential growth of 
the energy of the r-mode Hi 

In the above-mentioned investigations Hi , the lin¬ 
earized hydrodynamics equations with the GR force were 
used to obtain an expression for the time evolution of the 
physical energy of the r-mode perturbation, dE/dt, from 
which the gravitational radiation and viscous timescales 
were determined. In this work, in which the main results 
of Ref. 0 are reported, we present an explicit expression 
for the r-mode velocity perturbations that solves the lin¬ 
earized hydrodynamics equations with the GR force. Our 
conclusions are that (i) these velocity perturbations are 
sinusoidal with the same frequency as the well-known GR 
force-free linear r-mode solution, (ii) the GR force drives 
the r-modes unstable with a growth timescale that agrees 
with the expression found in Refs. Hi , and (iii) the 
amplitude of these velocity perturbations is corrected, 
relatively to the GR force-free case, by a term of order 
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dtv + {v-V)v = -p-^VP + (2) 

dtp + V ■{pv)=0, (3) 

V2$ = 47rGp, (4) 

where p, P and v are, respectively, the density, the pres¬ 
sure and the velocity of the fluid, $ is the Newtonian 
potential, and the GR reaction force, 

F^^ = -dt^+vxiVxp), (5) 

is assumed to be given by the 3.5 post-Newtonian order 
expansion that includes the contribution of the current 
quadrupole moment, which is the main responsible for 
the GR instability that sets in @. In Eq. 10 , ,5 is 
the gravitational vector potencial whose components are 
given by 


(^i = (6) 

where Sij (t) is the time-varying current quadrupole ten¬ 
sor, 


= J d^Xekq{iXj)XkpVq, (7) 

Cijk is the Levi-Givita tensor, Xi is the Gartesian coordi¬ 
nate of the point at which the tensor is evaluated, and 
denotes the time derivative of Sij. 
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III. LINEARIZED HYDRODYNAMICS 
EQUATIONS WITH GR REACTION FORCE 

The hydrodynamic equations ©-0 can be linearized, 
yielding 

dtS^^^v + ■ V)v + {v ■ V)S^^'>v 

= ( 8 ) 

^ • V(5(iV + V • = 0, (9) 

V25(i)$ = 47rG5(iV, ( 10 ) 


Using Eqs. m and 03, the first-order Eulerian 
change in the gravitational vector potential, is then 

computed to be 

= 0, (15a) 

SWpg = -K4aQa;^r2sin6»e-"^*e*(^'^+“‘>*), (15b) 

R 

s<~^% = -iK^afluj^r^ sine cos ee-^^F^‘^’^+‘^°^\{15c) 

R 

where the constant k sets the strength of the GR reaction 
force and is dehned as 


where v = Or sin 9e(j, is the velocity of the unperturbed 
star, p its mass density, denotes the first-order 

Eulerian change in a quantity Q and we have defined 

To compute the first-order Eulerian change in the GR 
force, (5 (i)f1^ii, we assume that the GR force-free r-mode 
velocity perturbationsi, 

(5(1^ = 0 , ( 11 a) 

= _l./I^^ 2 sin 0 e*( 2 <^+“‘)^ ( 11 b) 

4 V tt it 

sine cos ( 11 c) 

4 V tt it 

act as a source for the first-order Eulerian change in the 
current quadrupole tensor, 6'^^'iSij. In the above expres¬ 
sion, a is the amplitude of the r-mode and we assume 
that 


f^G 

45V 5 


(16) 


Finally, taking into account that Sij = 0, implying 
that 0i = 0, the first-order Eulerian change in the GR 
force, (5 (i)f1^ii, is computed to be 


j(1)^GR 


j(i)irGR 


5 ( 1 )^,GR 


—SjK-^crQ^WQr^ sin^ 9 cos 9 

R 

Xg-''^*gh 20 +‘^oi) 

iK^aVtujQ (wo -f 3Q sin^ 0) r^ sin 9 

^g-rotgi(2<^+(x,oi)^ 

—K^aflLOgT^ sin 9 cos 9 

R 

^g-rotgi(2<^+(x,ot) 


(17a) 


(17b) 


(17c) 


uj = ljq -\- ivj, (12) 

where the frequency of the r-mode, ojq = Re[w], and 
the small imaginary part that is related to the growth 
timescale of the instability of the mode, zu = Im[w] < 0 , 
are arbitrary parameters to be determined. 

Under the above assumptions, is computed to 

be 

<5(^(5,, = -aU./f 4 (13) 

\ 0 R 

where we have neglected the contribution coming from 
the terms ViS^^'^p (of order afl^) and retained only the 
dominant terms pS^^'^Vi (of order aU). Similarly, it is 
straightforward to show that the first-order Eulerian 
change in the other components of the quadrupole tensor 
satisfy the relations 

xWq — _a(i)q 

G ^yy — ^ ^xx 5 

= <5(1)5^, = 5(1)5,, = 0. (14) 


1 In this article, we will be concerned exclusively with the 1 = 2 
r-mode, which is the most unstable mode. 


IV. THE ANALYTICAL R-MODE SOLUTION 
WITH GR REACTION FORCE 


The linearized hydrodynamic equations (ISll-ltTUl. with 
the first-order Eulerian change in the GR force, 5 (i)f1^ii, 
given by Eq. ini), admit the solution 0 


5 ( 1 )?;^ 

= 0 , 


= -laU 


xe-“* 


= iaU 


xe~^* 


1 /5 1 


2 V 7ri? 


+ i^ [A — 1 ) fl® 




(18a) 
sin 9 
(18b) 
sin 9 cos 9 
(18c) 


and 


5 ( 1 ) 7 / = -aQ2 
3 


1 /5 1 


2 V 7ri? 


-I- 17 AQ® I sin^ 9 cos 9 


^g-CTtgi(2</)-|-aJot) 


(19) 


with zu and ojo given by 

8 /tF 


4fl 


= —A —jRi}° and ujq = -. ( 20 ) 

9 V 5 
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The velocity perturbations given by Eq. have a 
piece similar to the GR force-free solution, the differ¬ 
ence being the factor responsible for the exponen¬ 
tial growth of the r-mode amplitude due to the presence 
of a GR reaction force, and another piece proportional to 
a'y{A — l)n®, where A is a constant fixed by the choice 
of initial data and 7 = 1024/cJ/(81i?). 

The GR force-free limit is obtained when we set the 
parameter k, dehned in Eq. equal to zero. In this 

limit, w and 7 also go to zero. Then, from Eqs. ra and 
TO . we recover the GR force-free linear r-mode solution. 

V. DISCUSSION OF THE RESULTS AND 
FUTURE DIRECTIONS 

In this work we have presented the analytical r-mode 
solution to the linearized Newtonian hydrodynamic equa¬ 
tions with the GR reaction force. The velocity pertur¬ 
bations are proportional to with wq = 

—4f2/3. Thus, they have the same sinusoidal behavior 
and the same frequency loq as the GR force-free veloc¬ 
ity perturbations given by Eq. CD- The amplitude of 
the velocity perturbations is proportional to exp{—-ccrt}. 
Since tu < 0, the GR force induces then an exponential 
growth in the r-mode amplitude. The e-folding growth 
timescale tgr = l/w agrees with the GR timescale o 
found in Refs. The velocity perturbations 

contain also a piece proportional to 07 (A — 1)14®, where 
A is an arbitrary constant fixed by the choice of initial 
data. If we choose this constant A to be of order unity, 
then this part of the solution could be neglected 0 . 

A quite interesting feature that has emerged from re¬ 
cent investigations on r-modes is the presence of differen¬ 
tial rotation induced by the r-mode oscillation in a back¬ 


ground star that is initially uniformly rotating. That 
differential rotation drifts of kinematical nature could be 
induced by r-mode oscillations of the stellar fluid was first 
suggested in Ref. 0- The existence of these drifts was 
confirmed in numerical simulations of nonlinear r-modes 
carried out both in general relativistic hydrodynamics 
0 and in Newtonian hydrodynamics 0- Differential ro¬ 
tation was also reported in a model of a thin spherical 
shell of a rotating incompressible fluid M- Recently, an 
analytical solution, representing differential rotation of 
r-modes that produce large scale drifts of fluid elements 
along stellar latitudes, was found within the nonlinear 
Newtonian theory up to second order in the mode ampli¬ 
tude and in the absence of GR reaction M- This differ¬ 
ential rotation plays a relevant role in the nonlinear evo¬ 
lution of the r-mode instability 0 Two questions could 
be then naturally raised, namely, is differential rotation 
also induced by the gravitational radiation reaction and 
does this differential rotation play a relevant role in the 
nonlinear evolution of the r-mode instability? One of the 
aims of the investigation carried out in Ref. 0 is to initi¬ 
ate a programme that hopefully will allow to answer this 
question. The natural continuation of this investigation 
is then to try to find an analytical r-mode solution of the 
nonlinear hydrodynamic equations with the GR reaction 
force. This work is in progress 0 - 
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